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YAW STABILIZATION OF LANDING CRAFT 


RICHARD JON HINKLE 


Submitted to the Department of Naval Architecture and Marine 
Engineering on May 20, 1966 in partial fulfillment of the 
requirements for the Master of Science degree in Naval Archi- 
tecture and Marine Engineering. 

The purpose of this work was to investigate the applica- 
tion of pump jet propulsion and yaw angle stabilization for 
a landing craft going through the surf. Linear theory was used 
throughout the analysis of the control problem. The stability 
indices and derivatives were determined in part by experiment 
and the remainder by analytic techniques. 

It was found that LCVP landing craft are dynamically 
unstable without a skeg. Stability could possibly be restored 
by a feedback control system, but the far more practical solu- 
tion is to start with a dynamically stable craft. This will 
provide the craft with a greater margin of stability in all cases. 

Although it is felt that the craft can not be satisfactorily 
used without a skeg, much of the analytic data in Appendix A 
may be used to cover the case with a skeg. Reference (3) con- 
tains the appropriate modifications needed for the changes. 
Thesis Supervisor: MARTIN A. ABKOWITZ 


Title: Professor of Naval Architecture 
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velocity of hull cemter of gravity in the forward direction 
differential operator 
velocity of hull C.G. in starboard direction 


longitudinal distance from L.C.G. to the C.G. of the 
lateral added mass 


distance to L.C.G. from the reference point for forces, 
moments, velocities etc. 


Longitudinal distance from L.C.t%, to the center of 
pressure at which the lateral force Y¥ acts (taken as 
center of area of hull profiie) 

approximated by half the prismatic coefficient 

force in the starboard direction 

mass density of water 


yaw angle 


subscripts meaning partial derivative with respect to 
r&v respectively 


subscript meaning maximum value 
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INTRODUCTION 


Landing craft used in amphibious warfare are susceptable 
to broaching while transiting the surf. ‘The craft, im the surf, 
is ~articularly vulnerable because it i5 ina following sea 
and also the relative velocity of the water in the vicinity 
of the rudder is small, making the rudder less effective. 

A yossible method of reducing this broachiaog problem 
is to use a yaw angle control system. By replacing the screw 
proveller and rudder with a pump jet with a controlled dis-e 
charge angle, the problem of decreased rudder effectiveness 
would be eliminated. A pump jet ;ropulsion system of the 
Siz2 required for the craft weeds some meaus of power assist 
to vector the jet stream for steering. ‘This vectoring system 
might easily be converted to an automatic steering control 
system to minimize the yaw angle and eliminate breaching. 

The *conomics of war require that in order to be effective 
the craft must be cheap so the coutrol system tiust be of 
reiatively simple desigu and accuate for small peric?: of 
tir.e such as transiting a surf line. 

The equations of motion for the craft are a pair of 
linearized second order differeutial equations relating yaw 
an: sway to the yawing moments and the transverse force. 

The soefficients for the various terms im these equations 
weve determiaed by either experiment on a modei LCVF or the 
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aneiytical appreach of Jacobs (2). The current LCVPs are 
equipped with a skeg for improved dynamic stability and 
protection of the propeller. The tests run on the model at 
the towing tank were all conducted without this skeg. 

The analysis of the contrel preblem was done using 
limear control theory. The criteria of stability for the 
system was either that of Nyquist er by root locus. A model 
of the system for the analog computer was made, but due to the 


instability the results are not shown. 
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PROCEDURE 


From the general linearized equations of ship motion of 
Abkowitz (1) the equations for yaw and sway were obtained by 
decoupling the equations for the other modes of freedom. The 


resulting equations are: 


eee Re: 2) gigi Diy 
(x =mx,)OE + (y_-ma, ot, (Ym) + Fi) ae ~ Ya (1) 
dt dt 
acy ay dv 
(N -I_) + (N omxu J—— + (N -mx.)2- + Nv = =N .=N (2) 
a ar << o.. vy 6G V rud d 
b dt b dt b dt , 
el 22 23 24 
where Y and N stand for the forces or moments caused by 
rud rud 


the steering device.and the subscript d for the disturbance. 
Noting that any force can be resolved into a force and a moment 
at some other location, the force, due to deflecting the pump 
jet, will become the ol 


and Ne at the reference position 


| a 


for the equations. N nud is simply Y ay multiplied by the distance 


| 
k between its place of application and the reference point. A 
Similar treatment applies to disturbance forces and moments. 


xeducing the pair of simultaneous equations to one equation in Y 


and Y force we have: 
Wilh, obo 5-8), bdo) + (br, B04b. b) We ub, oe B, a + 
21°137711°23 iy el lgres oe alr os gtP 
(4 .Bp9rP 24) 2)8] = Fou g[CPzsrKb 58 + (D5,-kb,,)]+ ¥,[(b,,-k", 8 


(bo ,,7k! b 1) 
“ J 5 








Th. block diagram representation of the ship is shown in figure - 
The determination of the stability derivatives was done by 
experiment where the necessary equipment was available and the 
remaining derivatives evaluated analytically. The evaluation of 
NY was done by towing the model in the tank at a small yaw angle 
and measuring the moment. The moment was measured by construct= 
ing a heave rod with the lewer eight inches machined round and 
instrumenting the round portion with a full bridge of strain 
gauges. The full bridge array allowed the measurement si torque 
without any effect due to the bending moment in the bar. The 
Output of this dynamometer was averaged by an integrator to ob- 
tain the result. Measurement of + oe was accomplished at the same 
time by mounting the thrust block athwartships in the model te 
measure the force in the Y direction. Runs were made at three 
separate yaw angles at each of four different speeds. These 
measurements were made at the longitudinal center of gravity in 
order to have several terms drop out of the equations. The 


recaining stability derivatives were calculated using the method 


se. forth in Estimation of Stability Derivatives :.. 


Various Ship Forms, and Comparison With Experimental Results by 
Jacobs. This method employs a strip theory technique with 


lateral added mass coefficient along with Lamb’s coefficients 
of accession to inertia and hydrofoil.theory. 


The use of pump jetprepulsion introduces an additional 
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FIGURE I 


Block Diagram of Ship 
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fe .tor in both the e. and N. derivatives. Due to the resziivel’ 
high velocity of the water jet, a rotational momentum change force 
become significant. ‘The force derivative is equal to twice the mass 
of the water in the duct times the velocity of the water jet. The 
moment derivative is derived by applying this force half way down 
the longitudinal lemgth of the discharge duct. For purposes of 
Calculations a ten foot duct with @ six inch diameter was assumed. 
The thrust of the pump jet was determined by scaling up model 
resistance data obtained in the towing tank. an abterpt waa 

meade to gather resistance data for the model in «+ surf but the 
waves would not break for a long enough distance to gather 

any data. 

Because pump jet propulsion was employed, the skeg was not 
needed for protection of the propeller and was not on the model 
during the tests, The elimination of the skeg made the model 
dynamically unstable in straight line motion. 

The automatic control system was chosem to be one of few 
pects, The quantities and y are measured by a free gyro and a 
reuse gyro. It may be acceptable to have accurac, | ye 
periods of time such as when going through the surf. In this case 
the measurement of # can be accomplished by using the gyro as a 
free gyroscepe rather than as a null device. If the gyroscope is 
frea, some means of sensing its positiom without applying much 
torque, such as a differential transformer, will have te be used to 
avoid unwanted precession. The values of ¥ and y times their appro= 


priate gains are sunmed and compared with the rudder angle and 
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the difference applied to a push pull amplifier. The output of 
the amplifier is a solenoid shown in figure Il. By Keeping the 
Magnetic field in the solenoid large enough the plunger will 
become magnetically saturated. The force exerted by the coils 
on the plunger will then become directly proportional to the 
current. The position of the plunger is used to control an under= 
lapped spool valve of the constant pressure hydraulic system. 
The reaction of the spool valve as a epring and damper were con- 
Sidered so that this portion of the system could be critically 
damped. This was done to help reduce the high fresuency oscilla- 
tions of the vectoring duct. The output from this four way valve 
drives the rams that position the jet vectoring duct. 

The removal of the skeg made the craft dynamically unstable, 
The next step was to investigate what could be done by the con- 
trol system to restore the dynamical stability. First, the effects 
of control, proportional only to yaw angle, was investigated. 
It was shown by the Nyquist Stability Criteria that, at the speed 
chosen for investigation, no amount of gain would suffice for 
Stabilization. Proceeding to a control proportional to yaw angle 
and yaw rate again, it could be shown that the system would be 
ineffective. Using a root locus plot, it can be seen that the 
right most two poles would continue to stay on the right half plane 
for all values of gain. By studing the root locus further, it 
can be seen that a system proportional to yaw plus the first and 


second derivatives would be the minimum requirement to produce 
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a dynamically stable craft by adding two complex comiuga’r cers 
to the root locus through the use of appropriate feed back gains, 
the unstable roots could be brought over to the ieft half plane. 
The design of such a system was dismissed however. nydill (3) 
claims that such systems are not satisfactory due to excessive 
high frequency oscillations of the steering device. Also, the 
expense of the extra equipment needed for the contrel geuld be 
eliminated by reinstalling the skeg and making the craft dynam- 


ically stable to begin with. 
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RESCLT 


The results of this investigation were based on the following 


assumptions regarding the design of the system: 


Ge 
De 
Ce 
d. 
Go 
f. 
ge 
he 


i. 


je 


longitudinal length of jet duct is 10 feet 
jet duct diameter is 6 inches 
Y force from deflected jet acts 15.94 feet from L.U.G. 


density of the hydraulic fluid 57 lbs/ft? 


2 
8=-SEC 
in 


mass of spool valve and plunger is 6.88 x io 
discharge coefficients for spool valve ~ : 
width of spool valve ports are 0.75 iiches 

the rams that drive the jet vectoring conduit are 1 
square inch in area ana are located 6 inches from the 
swivel, therefore the discharge angle is equal to the 
ram travel in inches over 6 


supply pressure is constant at 1000 psi 


the radius of gyration equals one forth the length 


The thrust to be developed must equal the resistance of the 


craft for constant velocity operation. The model was tested in the 


towing tank and the model resistance data and scaled up resi. 


are tabulated in figure III. A plot of the full scale resistance 


is shown in figure IV. 


The qualities Y and NY were found experimentally in the 


towing tank. The three yaw angles used during the testing were 


Dete®, 5.92°, and 9.30° and the model speeds varied from 1.3 knots 
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to 2.9 knots. ‘This corresponded to speeds between 4.75 and 10.5 
knots for the full scale craft. a table of the experimental data 
is presented in figure V. The data was non-dementionalized by 
dividing forces by hy u“1° and moments by [Pee Both NY and oe 
were found to be dependent on Froude Number as is shown by the 
plot in figure VI. ‘lhe variation of yaw angle due to the twist- 
ing of the rod was computed for the largest torque experienced 
and the effect was found to be insignificant, 

The data for the remainder of the stability derivatives and 

indices was computed by the method outlined in reference e,. 
The method of nonedementionalizing used here was to divide forces 
by W- The length used in all calculations here was that of the 
waterline 32,367 feet. ‘the actual computations of these deriva- 
tives and indices are contained in appendix A. 

The value of stability indices, derivatives, and resistance 
for speeds of 4 and 7.5 knots are listed in figure V1l. 

The spool valve is an underlapped four way valves with 
ports 0.75 inches wide and discharge coefficients of 0.7 at each 
port. The equations for flow rate for such a valve is that 
ja R\Pos where j is ram rate, X is spool value displacement, 

P is supply pressure, k is the product of discharge coefficient, 
port width and the square reot of two times the acceleration of 
gravity divided by the mass of the density of the hydraulic fluic, 
(xaew/ 2g- This equation gives j=3600x for this system. The 


spring reaction force of the valwe is given by the equation: 


15 
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a 


k=2 x discharge coefficient x port width x Pi x cos 09° 


g 


or k=2 x 0.7 x 0.75 x 10° x 0.358 = 376 ie 


The damping reaction of the spool valve is expressed by the 
equation: 

C=discharge coefficient x port width x Ves x L 
where L is the difference between distance from supply port to 
load port and load port to drain port. With L equal to 2, the 
Welwe of C is 0.7 x 0.75 x 05293 x 220.5380 


By using a spring with a constant of 24 ae to center the 


g lbs=sec 
in 


plunger, the total spring constant becomes 400 iS The plunger=_ 


2 
spool valve combination is assumed to have a mass of 6.88 x 16” Eieeeee 


n 
To critically damp this system, a damping constant of 3 
is needed. An additional damper with a constant of Ele 


would be needed to critically damp the systen. 
The combined control system and ship block diagram for 


speeds of 4 and 7.5 knots is illustrated in figures VIII. 
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FIGURE IIl 


RESISTANCE DATA 


Ne poy C, 
4.33. x 107 0.05 0.00974 
7.91 x 107 0.1575 0.00914 
1,061) & 10 0.283 0.00924 
15428 x 10° 0.668 0.01022 
1.750 x 10 1.293 0.01533 
2.115 x 10’ 2.413 0.01960 
2.282 x 10° 3.593 0.02515 

N. Ce +DC, C, 
8.07 x 10° 3.439 x 10° 7.97 x 
1.43 x 10° 3.170 x 10° 7.69 x 
1.97 x 10° 3,035 x 10° 7.92 x 


2.67 x 10° 2.913 x 10° 8.91 x 


Ce 


5.209 x 10° 4.53 x 16° 
4.615 x 10° 4.52 x 10° 
4.363 x 10° 4.88 x 10° 
4.119 x 107 6.00 x 10° 


20060 cuir Iidiels 107 


4,831 x 10? 1.577 x 10° 


3.779 x 10° 2.137 x 10° 


Resisteace 
los 


16° 21.25 
16° 69.1 
10? 126.4 


10° 260.5 


3.275 x 102.837 x 10° 1.400 x 10° 616 


3.94 x 10° 2.771 x 10° 1.854 x 10° 1182 


4.26 x 10" 2.743 x 10° 2.411 x 10° 1807 
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FIGURE V 


Experimental Values of Yyand Nv 








Y=3,12° 
i. tog Nin lbs = ¥? Ne 
2.27 0.0613 203 6.37 x 120° 4,32 x 10° 
2.27 0.0580 2.03 5.88x10' 3.81 x 10° 
3.03 0.0833 3.48 4.72 x10' 3.67 x 10° 
47? 004868 5015 5043 x 10° 26875 x 10° 
W=5.92° 
2.27 0.1162 3.70 1.173 x 10° 6.95 x 10° 
Ba@t =s—s«C« 1860S 5018 6680 x 15° Samp & 16¢ 
BeOS) 86=—s«a 1085755 BK x 16" 5.00 x 16¢ 
4o1? —0K66 BK 14625 x 10? 4,68 x 10° 
4.54 095307 1040 14342 x 10° 4.89 x 10° 
W9.30° 
2.055 Oslhl4 2.73 16749 x 10° 6.26 x 10° 
2.17 0.1818 3.75 1.83 x 10° 7.04 x 15" 
3203 0.3404 6.30 1.93 x10? 6.64 x 10° 
4417 047405 12010 2.22 x 10° 6.74 x 10° 
5.03 1.2020 20.40 2.47 x 10? 7.80 x 10° 
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FIGURE VII 


STABILITIES INDICES AND DERIVATIVES 


Coriolus Force 
Coriolus Moment 


Resistance 
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FIGURE VIII 


Block Diagram of Ship and Control System 


15° knots 
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DISCUSSION. OF RESULTS 


The instability of the system with control proportional only 
to yaw angle is shown in figure 1X, a Nyquist Stability Criteria 
plot at 7.5 knots. There is one pole of the open loop transfer 
function in the right half plane. For ynity. gain in the control 
system, there are no tiircieuenta of the minus one point, 
therefore there will be one root in the right half plane of the 
closed loop transfer fumction. B,; increasing the gain, the small 
loop near the origin (shown in the inset of figure IX) will 
expand and eventually encircle the minus one point once. This 
indicates that increasing the gain will only add to the instabil- 
ity by putting another pole in the right half plane. The Nyquist 
plot at a speed of 4 knots is not significantly different to 
change the results regarding stability. 

To add to the sophistication of the control system by includ-= 
yaw rate feedback in the system, we can most essily visualize the 
Stability problem by a root lecus plot. The determination of the 
exact breakaway point requires the solution of a sixth degree 
equetion, but the trend of the system can be seen without the 
exact solution. The poles aiid zeros of the open loop transfer 
function and the asymptotes for speed of 4 and 7.5 knots are 
plotted in figures X and XJ]. The effect of preportional plus 
first derivative control is to add another zero to the plot and 
1ts poSition would be governed by the ratio of the gains of the 


proportional te the Sirst derivative terms. With the total 
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system gain at zero, we will already have one pole of the closed 
loop transfer in the right half plane. As the gain is increased, 
the rightmost two poles will come together, merge, and split, 
injicating complex conjugate roots. The roots would follow wene® 
asymptotes originating from approximately the -95.$ point. At 
no time would all the roots be in the left half plane. The placing 
of a zero at the origin by this method of feedback control, could 
only be accomplished with infinite gaiu on the derivative feed- 
back and so is dismissed as a solution. The only alternative will 
be the speeding up of the coutrol system, such ae by reducing 
the mass of the plunger and spool valve. This would in effect 
move the first pole on the left of the origin further left. 
After the pole passes to the left of the zero, it may become 
possible, with sufficiently high gains, to bring all the roots 
to the left half plane and stabilize the system. 

Proceeding to the next step would be to make the control 
Sy«tem proportional plus first and second derivative control. 
Hert again, the simplest way to visualize the stability problem 
is with a root locus plet. ‘lhe effect of such a svete” would be 
to add two zeros to the plet. By selection of appropriate gain 
ra.i0s, these zeros will cause the poles on the right half plane 
to cross the imaginary axis at a relatively low value of gain. 
Th. asymptotes in this case will be at 60° above and below the 
re . axis and on the negative real axis. The point of inter- 
section of the asymptetes and the real axis will be governed by the 
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yeal par’ o1 the position of the two new zeros. The 60° slopes 
obviously indicates that excessive gain will again place the roots 
in the right half plane and unstabilize the system. As was 
mentioned before, systems with second derivative feedback have 
not proved satisfactory for automatic steering control aboard 
ship. ihe control system tends to produce relatively high 
frequency oscillations of the steering device. A continual meve= 
ment of the @teearing control would be necessary to stabilize a 
dynumiraliy ulisteble ship, but these systems tend to cause much 
more wotinm thang 1s needed fer effective control. The introduce 
tion of a second degree lag network into the control system is 


fediiweeedead bey w¥dat! (4) ta wel: reduce these oscillations. 








CONCLUSION 


.t 1s evident that any attempt at automatic control for yaw 

stub lization of an LOVE without a skeg, will ke difficult and 
uneconomical. By using the skeg to restore dynamic stability 

to the craft before the control system is added, a great deal of 
Coupderaity will be eliminated, The dynamically stable craft will 
naturally be ®asier to handle wher steering manually. The use of 
1 es ion appesls to be a practical method of prepelling 
svt contreliling landing craft and is adaptable to automatic 


cantrol, 


oe 











RECOMMENDATIONS 


‘The most importamt recommendation can only be to use a skeg 
to make the LOV: dynamically stable to begin with. If this 
scheme of testing is to be used again for this model, Jacob's 
method may be used to include the use of the skeg too. 

The tewting in the towing tank should include some means 

f obtaining tesistance data while in a surf. This will govern 

fe thrust avo1rlanle for control. cs suggest that some ramp be 
installLe’ in the tanh near the wave maker to generate breakers. 
shio ramp «ould have the same effect as a reef in causing the 
wavwee tu tecome breakérs. IJ feel any attempts to gather data 


iremkiog im the tank's beach will prove futile, 


due tw the short length of surf, 


28 











BIBLIOGRAPHY 


1. Abkowitz, M. A., Lectures 2n Ship Hydrodynamics-Steering 


and Mameuverability, Danish Technical Press, May 1964. 


2. vJacobs, We. Ro, Estimation of Stability Derivatives and 
Indices of Various Ship Forms, and Comparison with Ex=- 
perimental Results, Davidson Laboratory Report 1035, 
September, 1964, 


5e Rydill, L. J., A Linear Theory for the Steered Motion of 
Ships in Jdaves, Transaction of the Institution of Naval 
Architects, Vol. 101, 1959, pp. 81-112. 


4, Sehiff, L. J. and Gimpirch, M., sutomatic Steering of 
Ships by Froportional Control, Sogiety of Naval Archi- 
tects and Marine Engineers Transactions, Vol. 57, 1949, 
pp. 942124, 


29 























AFERNDIX A 


The following is a prref description of the Jacobs’ (2) method 

of determining stability indices. The first step was to deter- 
mine the sectional inertia coefficients C by entering figure XVII 
with the values cf sectional area divided by local beam times 

a ae ; * : | b~ 
loosed drwrt Me? and luceal beam divided by local draft “ 
Neat, inab's cwoeffariernts of accemsion to inertia (kek, ok*) 

2 5 < 4 r ops on ° Q 
ace found from figure dAvial by entering wath Le The prime in 
the following squaticws will denote aon-dementionalized quantities. 
(\GOte ali juantities here 2re nou-dementionalized by the quanti- 


ties noted is: tne charter on results) The limits on the 


INLEBLALS A+R cetressius buw to stern, 


ra 
? AD -, 
® 2h. 4 2 Gb da” (4) 


Y= ~<a 
e. An (5) 


r 
> it a 
Ly: 4 (6) 


NI 





§ 
atm! 
ak, 2 (7) 
re 
j, Cb Base 
i ~_ 
ke 2 (3) 
—% Ch de 
cD 
% ‘ ’ or: 
~aphe: nal{ the [risomatacs coefficient (9) 


od 








yttF 22 
Ne =k! PAM o on x dx (10) 
Ye sm) =k, Mm, (11) 
Yo =-m,X (12) 
Ne ==m_x (18) 
2 
f ’ 0 
Yi skint + ms )-5 Ly (14) 
iyi @ 
Y =, /anl (15) 
N° ee (16) 
r ahi a 
=nt Wuni® 
Np =Nt Aunl (17) 


Figures XII through XVI are the calculations required to 


derive these stability derivatives and indices for the LCVP 





tested. The integrations were carried out by Simpson's rule 


and all areas found by use of a planimeter. 
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FIGURE AIT 
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2 0 O 0 O 4 0 O 
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ed 3.50 0.97 3.40 +15.09 1 3.40 =51.3 
134.73 38.17 

¥= Pc n®xax = gli = 0.2835 

uA = ix Bh 
Mp 
han, '/, a dX2G 0985 x 2252 Ps 134.73=39305 





FIGURE XII (continued) 


ter : 
Ye=k 5 ra Ch xdx=0.935 x 0.995 x 3.14159 x 38.17=111.5 











FIGURE XITI 
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FIGURE aALV 
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FIGURE XV 
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FIGURE XVI 


Ns =0.813 x 2202 , 7443 = 18,890 


T=7e14 x 65.4 x 10°s46,700 
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